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As an example of ordering due to quantum fluctuations, we examine the nearest-neighbor antifer- 
romagnetic quantum 0{n) rotor model on the pyrochlore lattice. Classically, this system remains 
disordered even at zero temperature; we find that adding quantum fluctuations induces an ordered 
phase that survives to positive temperature, and we determine how its phase diagram scales with 
the coupling constant and the number of spin components. We demonstrate, using quantum Monte 
Carlo simulations, that this phase has long-range spin-nematic order, and that the phase transition 
into it appears to be flrst order. 
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I. INTRODUCTION 

Ordering of antiferromagnetic spins on geometrically 
frustrated lattices is a subtle problem P, • By defi- 
nition, the canonical nearest neighbor antiferromagnetic 
interaction on such lattices fails to produce a unique clas- 
sical ground state. Thus low temperature ordering in 
such systems must be attributed to additional interac- 
tions or to selection via thermal or quantum fluctuations, 
the latter phenomenon being termed order by disorder. 

In recent years there has been much interest in anti- 
ferromagnetism on the pyrochlore lattice. Figure 
The classical statistical mechanics of the purely nearest 
neighbor problem is now fairly well understood. It is 
known 's'l that XY spins order, by disorder, coUinearly 
[SiJ while Ising, Heisenberg and higher-dimensional spins 
do not order even in the zero temperature limit. The 
correlations of this set of cooperative paramagnets have 
been found to exhibit a universal dipolar form charac- 
teristic of an underlying gauge field 6;|, truncated by a 
correlation length that diverges as T ^ 0. 

While thermal fluctuations thus do not lead to order- 
ing in the Heisenberg problem, there is much work ar- 
guing that quantum fluctuations do lead to ordering at 
low temperatures, for the case of Heisenberg spins. Close 
to the classical limit, this conclusion follows from argu- 
ments based on the 1/S expansion discussed recently by 
Henley 0- Henley derives an effective Hamiltonian on 
the space of classical ground states which captures the 
effects of the zero point energy of harmonic spin waves in 
a loop expansion. This indicates a selection of coUinear 
ground states with a large unit cell and a residual de- 
generacy that is of order 0{L), where L is the linear size 
of the sample. It is expected that a nonlinear treatment 
of the spinwaves will lift this remaining degeneracy and 
predict long range spin order in a particular coUinear 
configuration. In the opposite limit of small spins and 
large quantum fluctations, there is a set of investigations 
principally of the 5=1/2 case, starting with the pio- 
neering work of Harris et al. [a, 0, ITol Ull which 
suggest a breaking of the inversion and translation sym- 
metries of the lattice with only short range order among 
the spins. For 5 = 1 Yamashita et al.^^ have also ar- 



gued for a breaking of inversion symmetry as well as long 
range order in the transverse component of the spin chi- 
rality. Two large n studies 15, do] have also suggested 
symmetry breaking for small spin values. Indeed Ref. [T^ 
finds that the quantum dimer model that arises in a large 
N treatment, does indeed break inversion symmetry and 
further translational symmetry breaking via an order by 
disorder mechanism. 

In this paper we further explore the impact of quan- 
tum fluctuations on spins on the pyrochlore lattice by 
endowing them with the dynamics of quantum rotors in- 
stead. The resulting problems, which are readily defined 
for 0{n) symmetric spins for all n, are distinct from the 
Heisenberg spin problems even at n = 3 although there is 
clearly a family resemblance which detailed analysis will 
bear out. The rotor models are interesting in their own 
right and their n — 1 representative is the transverse 
field Ising model which has been studied on a variety of 
frustrated lattices [T7| . 




FIG. 1: Pyrochlore Lattice 

The Hamiltonian for the quantum rotor model on the 
pyrochlore lattice is given by 



H 



(1) 



<«j> 



where Si is a unit vector in n dimensions located at site 
i, Li is its associated angular momentum and the second 
sum goes over all nearest-neighbor pairs of sites on the 
pyrochlore lattice. The coupling constant g measures the 
strength of quantum fluctuations. 
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At any n and in the limit of large coupling constan 
g, the ground state has all rotors in their zero angula 
momentum state, and thus no long range correlations o 
symmetry breaking. In this limit there is also a gap o 
0{g) so this quantum paramagnetic phase is stable fo 
a range of values of g. [IJ The opposite limit of g = 

is the classical model, which has a highly degenerati 
ground state and, for n > 2, dipolar correlations in thi 
zero-temperature limit but no significant strength in am 
Fourier component of the magnetization |^ |^ . 

The problem at hand is to understand how the sys 
tern interpolates between these two very different disor 
dered limits. In the following we report progress on thi; 
question. Primarily we will show that the system devel 
ops spin-nematic (collinear) long-range order for a finiti 
range of values of g e (0, gc) at T = 0. More generally wi 
derive the schematic phase diagram in the temperature 
coupling constant plane shown in Fig. |21 with the var 
ious scalings that we have derived indicated. We hav 
not been able to establish or exclude further symmetry 
breaking at the lowest temperatures into a state with long 
range spin order. We do however offer analytic evidence 
and simulational evidence (for the case n = 3) that such 
further symmetry breaking is likely quite weak. 

We turn next to understanding the quantum mechan- 
ics of the building block of the pyrochlore lattice, the 
single tetrahedron. The basic results derived here will 
next enable us to deduce the existence of nematic order 
on the full lattice and the scalings of various energy scales 
with n. Thereafter we describe results of a simulation for 
the n — 3 problem and conclude with a discussion of the 

1 /n expansion for this problem and a summary. 

II. SINGLE TETRAHEDRON 

In this section we analyze the basic unit of the py- 
rochlore lattice — a single tetrahedron. It consists of 
four mutually coupled spins and its hamiltonian is 

4 

H = gY,Lj+ ^^-^i- (2) 

i=l l<i<i<4 

This may also be viewed as a system of four interacting 
particles, each moving on the surface of a sphere in n 
dimensions. The first term in H is their kinetic energy, 
while the second term is a repulsive interactionpotential. 

In the classical limit of g = it is known [3 that for 
n = 2 (XY spins) the spins order coUinearly in the zero T 
limit — two spins pointing in one direction and the other 
two in the opposite direction On the other hand, for 
n > 2 the spins remain disordered even in the T ^ limit 
as we review below. In the rest of this section we show 
that in the quantum case {g > 0) for any n there are ways 
to take the g ^ and T ^ limits that give collinearly 
ordered spins. Furthermore, we find the dependence of 
the crossover from ordered to disordered spin states as a 
function of g, T and n, as indicated in Figure|21 
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FIG. 2: Schematic crossover diagram of a single tetrahedron, 
with the scalings with g, T and n indicated. The solid line is 
where the system crosses over from being collinear to disor- 
dered, while the dashed line shows the value of g where, for a 
given T, the deviations from coUinearity are minimized. The 
solid line is also a very approximate phase transition curve for 
the whole system. However we expect that the curve starting 
at gc will bend initially to the right due to the entropy of the 
spin waves that exist in the ordered phase. 



A. Classical spins, g = 

We start by describing the configurations of the spins 
and the classical ground states (cgs). The spins are in 
a cgs if the potential is minimized; this is when the four 
spins add up to zero. This implies that in a cgs, the 
fourth spin lies in the three-dimensional space spanned 
by the other three spins. Thus the ground states can be 
parameterized by two angles Q and as shown in Figure 
|31 plus an overall rotation in spin space. Q gives the 
deviation of all spins from coUinearity, while is the angle 
between the planes in spin space spanned by each pair of 
spins. All four spins in a cgs are at the same angle d from 
the reference axis; the reference axis that minimizes Q is 
used in our convention of parameterizing these ground 
states. Mostly we will consider nearly collinear states 
with small d. 

For both the classical and the quantum analysis, we 
need to consider the matrix of second derivatives of the 
potential energy with respect to the orientations of the 
four spins at a given classical ground state. First, this 
matrix generically has (3n — 4) eigenmodes with zero sec- 
ond derivative; these directions span the local manifold 
of classical ground states. For small Q it also has (n — 1) 
eigenmodes with second derivative of order unity, and 
one "soft mode" with second derivative proportional to 
(fl . This last mode corresponds to the spin deviation 
away from the nearly collinear cgs that decreases Q for 
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FIG. 3: a) Parametrization of the ground state, showing th' 
four spins and the reference axis. All the spins are at angle ( 
from the reference axis. The short arrows indicate the modi 
that goes soft as the spins become coUinear. The line with n( 
arrows is the reference axis in spin space that we measure ( 
from. 

one pair of nearly parallel spins and increases it for thi 
other pair, as indicated in Fig. 2. This soft mode's stiff 
ness is independent of the other angle (p. 

At low temperature, the classical system at equilibriun 
can be approximated as occupying the states with potcn 
tial energy within fcsT of the ground state and not thosi 
at higher potential. Due to the soft mode, the number o 
such states near an almost coUinear cgs is proportional 
to T/e (for T kO'^). For n = 2, this concentration of the 
accessible states near = 0, due to the one mode that 
softens there, causes the "order by disorder" effect and 
the spins order coUinearly in the T — limit. But the 
number of cgs with a given is proportional to 02(n-2)^ 
due to the freedom of rotations about the reference axis. 
This means that the probability density of behaves as 
~ 5((2n-5) g^i- sjj^all 0. Thus for n > 2 the coUinear states 
do not dominate even in the zero temperature limit, and 
the classical spins remain disordered. 



B. Quantum ground state, g > 

Here we show that for small g the ground state wave- 
function is localized around the coUinear state and we 
obtain the scaling with g and n of its spread, its energy 
and the energy of the lowest excited states. First we 
present a "power-counting" variational argument 

We treat the quantum zero-point energy of the mo- 
tion normal to the manifold of classical ground states as 
that of harmonic oscillators. For small 0, there are (n— 1) 
"stiff" oscillators, each with zero point energy ~ and 
the one soft mode, with zero point energy ~ 0^/g- The 
latter produces a 0-dependent effective potential that lifts 
the degeneracy within the set of classical ground states 




FIG. 4: Very schematic drawing of evolution of the distribu- 
tion of the eigenstates and probability of the system at a fixed 
T as <; is increased. The height of each peak corresponds to 
the probability that the system is at the given 6. The width of 
each peak corresponds to the spread of each eigenstate. The 
distribution is localized approximately up to an angle Of and 
we say that the system is localized if 0/ < 0; ~ 1. At the top 
g is small, the system is noncoUinear and effectively classical, 
occupying many excited states. As g is increased (moving 
down in the figure) the system becomes more coUinear, and 
the number of different excited states occupied at equilibrium 
decreases. Meanwhile, the spread of 6 within the ground state 
is always increasing as g increases. At pmin the excitation 
energy of the lowest excited state passes through T; this is 
where the system is most coUinear. As g increases further, 
the spread within the ground state increases (bottom). 



and has a minimum aX 9 — Q. We are interested in char- 
acterizing the ground state and the excited states in this 
potential. 

The set of classical ground states for a given reference 
axis constitutes a (2n — 3)-dimensional manifold that in- 
tersects with itself at the coUinear states where 6* = 0. 
The angle 6 is proportional to the distance from this in- 
tersection. Each pair of nearly parallel spins may be ro- 
tated about the reference axis. These latter rotations give 
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2(n — 2) dimensions of motion. One can also rotate the 
reference axis, which gives the remaining (n — 1) dimen- 
sions of motion. In the ground state, these reference- axis 
degrees of freedom are in the zero total angular momen- 
tum eigenstate and we will generally ignore them in this 
section of the paper. 

Thus we consider a variational wavcfunction that is 
localized near = 0, spreading in angle by 7 in each di- 
rection on the (2n — 3)-dimcnsional manifolds near the 
collinear state. The kinetic energy of such a state is 
^ g{2n — 3)/7^. The typical angle in this state is the 
(Pythagorean) sum of (n— 1) angles of displacement away 
from the reference axis that are mutually perpendicular 
to each other and each of order 7, so ^ ~ jy/n — 1. At the 
level of accuracy we are using now, n ~ (n — 1) ~ (2n— 3), 
so the kinetic energy is ~ and the full variational 

energy is this plus the effective potential of ^ Q^fg- Mini- 
mizing this with respect to B gives the following estimate 
for its typical value in the ground state: 



(3) 



Thus we see that the ground state is collinear (^0 0) 
in the 5 — > limit. The crossover from collinear ordering 

(small ^o) to a strongly disordered state can be defined 
as occurring at ~ 1) which puts the crossover in the 
ground state at 



n 



(4) 



for large n. 

The contributions to the ground state energy from the 
motion within the manifold of classical ground states as 
well as that from motion along the "soft mode" direction 
are 



(ngp 



(5) 



In the former case, this energy is due to the motion in 
(2n — 3) directions, so the energy of motion in just one 
of these directions along the manifolds of cgs, and thus 
the energy to excite the motion in that direction to a 
higher-energy eigenstate is 



(6) 



There are other, lower-lying excited states that do not al- 
ter the degree of collinearity: these involve "rigid body" 
rotations of the ground state we are discussing, and, at 
even lower energy for small g, tunneling between the 
three distinct ways of pairing the four spins. 

The above arguments have approximated the soft 
mode as harmonic, which is correct for the bulk of the 
ground state wavefunction in the limit of large n. How- 
ever, we should check that the g-dependence that we have 
derived remains valid at small g even when n is small 
(such as for the interesting case of n = 3). To do bet- 
ter at small g and small 6 but general n, we combine 
the soft mode degree of freedom with the manifold of 
classical ground states, and formulate the Schrodinger 



equation within this space. Label the two pairs of nearly 
parallel spins 1 and 2 and let the displacement of one 
of the spins in each pair away from the reference axis 
be fi and f'2, respectively. The other spin in each pair 
has precisely the opposite displacement. These displace- 
ments are (n — l)-dimensional vectors and, when small, 
their magnitudes are equal to the angles of rotation of 
the unit vectors. To leading order at small g and small 
ri the resulting Schrodinger equation is 

H'i, = -2g{Vl + Vi)V + (l/2)(r? - ri)V = E'^, (7) 

By comparing the kinetic and potential energy terms 
in this Hamiltonian H', we see that indeed the charac- 
teristic scale of angle is 56 and the scale of energy is 
53, as obtained above in the harmonic approximation. 
The ground state in this unusual quartic potential for 
large n is concentrated away from the collinear state at 
n = r2 ^ ns ge , where the harmonic approximation used 
before remains valid. But for small n the ground state 
wavefunction has considerable weight close to collinear- 
ity where the harmonic approximation is not appropri- 
ate. The leading correction at small g to (7) appears to be 
from the ^-dependence of the stiffness of the other "hard" 
modes that were ignored. This gives a ^ \fg contribu- 
tion to the effective potential that thus contributes to 
the energy at order g"^ , one order higher in our "small" 
parameter of g(6 . 



C. T>0 

Next we consider nonzero temperature for small 5, ex- 
amining the crossovers that occur, first from the fully 
quantum regime for T < E^,j. where the system remains 
in its nearly-coUinear ground state, to an intermediate 
regime (£Jex < T < Tc{g) ~ g^ /{n — 2)) where there 
are many thermal excitations present but it remains near 
collinear, and then, for n > 2, to the disordered regime 
at higher T. 

The modes within the manifold of classical ground 

states are excited when the temperature reaches the ex- 

2 1 

citation energy Tex = E^x g^ /n^ . At higher tempera- 
tures, we can treat these degrees of freedom as classical. 
The "soft mode" has an excitation energy of ^ 6^ and 
this mode will be in its classical regime where it is highly 
excited only at angles where this is less than T. There 
the probability of being near a particular cgs is 



P{0) 



0^9 ' 



(8) 



this applies for 9q < 9 < Tj For 9 < 9a, within the 
support of the ground state wavefunction, P{9) ^ P{9q). 
At larger angles, greater than both 9q and T/ the 
soft mode is in its ground state, which gives an effective 
energy of ~ 9^ and 



P{e)^e^p{-0^/T) . 



(9) 
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Again, the number of distinct cgs with a given be- 
haves as ^ Q'^i'^-'^) at small 6, so the typical value of 6 is 
near the maximum of 6l2("-2)p(-0) ^jjg intermediate 
regime we are discussing, which is E^x < T < Tdg) ~ 
f/2 /(n — 2), this maximum occurs near 

OiT) ^ . (10) 

Thus, for 71 > 2 and small T and g, the crossover to the 
thermally disordered state occurs at 

T,{g)^g-^/{n-2) , (11) 

where 0{T) ^ 1. These small g, T results apply as long 
as the quantum ground state is itself nearly coUinear, 
which requires g < and thus T < n~^. The crossover 
temperature Tc{g) must have a maximum of order 
and then decrease to zero at gc, as indicated in Fig. 1. 

The striking result here is that this simple system of 
four rotors has a nonmonotonic, or re-entrant behavior 
at low T as one increases g from the classical limit of 
g = to the quantum limit of large g for n > 2. This 
is illustrated in Fig. 3. Initially at small g it is disor- 
dered due to the large entropy of the disordered states 
relative to the coUinear states. As g is increased, the ef- 
fective potential due to the soft mode increases, causing 
the system to be more and more confined to the nearly 
coUinear eigenstates, as 0{T) decreases with increasing g. 
This trend continues until the energy E^x of the excited 
states within the manifold of cgs increases to of order T, 
at which point the system is predominantly in the ground 
state, with 0o ~ Tin^. This point, where the deviations 
from coUinearity are minimized, occurs at f/min ~ T^n^ . 
Beyond this point, further increase of g (decrease of the 
"mass" ) causes the ground state to instead deviate more 
from coUinearity with increasing g^ until it crosses over 
in to the fully disordered quantum regime at gc- 

III. PYROCHLORE LATTICE 

Now that we have examined the behavior of a single 
tetrahedron, we turn to the question of the behavior of 
our model on the pyrochlore lattice. This lattice consists 
of a three-dimensional array of corner-sharing tetrahe- 
dra, so adjacent tetrahedra share a single site. There 
are two aspects of the ordering that occurs in the single 
tetrahedron in the appropriately-taken small T limit, 
namely the axis along which the spins are all coUinearly 
aligned, and which pairs of spins are pointing which way 
along that axis. If two adjacent tetrahedra sharing a sin- 
gle spin both order, they must order along the same axis, 
so this "spin-nematic" order should propagate through- 
out the lattice. As we show below, we have found good 
evidence from quantum Monte Carlo simulations that 
the pyrochlore lattice model has a phase with long-range 
spin-nematic order in a region of its phase diagram with 
nonzero g and T . Thus the crossovers we discussed in 



the previous section for the single tetrahedron become 
true phase transitions on the full lattice. Since this is 
an isotropic-to-nematic phase transition with a cubic in- 
variant in its Landau theory, this phase transition is ex- 
pected to be first order, and that expectation is indeed 
supported by our simulations. 

The other aspect of the ordering does not strongly 
propagate between adjacent tetrahedra: The spin ne- 
matic order picks a particular axis in spin space, but 
there are 6 different low-energy configurations of the 
spins' directions along this axis for each tetrahedron. If a 
given tetrahedron orders in to one of these spin patterns, 
the adjacent tetrahedra that each share a spin with it 
still each have 3 spin patterns that remain compatible 
with the first tetrahedron. At this level of consideration, 
for the full lattice the entropy of this spin degeneracy 
is extensive, and the spins remain disordered, although 
coUinear. It seems likely that this spin degeneracy is 
lifted to some degree by tunneling between the various 
spin configurations, but this is something that so far we 
have not detected, either analytically or in our simula- 
tions. 



A. Degrees of freedom and modes 

We start by reviewing the counting of the number of 
degrees of freedom and of constraints j^. Let N be the 
number of tetrahedra. Then there are 2N spins and 
2N{n — 1) degrees of freedom. Each tetrahedron has zero 
total spin in a classical ground state; this gives Nn con- 
straints. Thus the dimension of the set of cgs is N{n—2), 
and hence there are N{n—2) zero modes around a generic 
cgs configuration. However, around a fully coUinear con- 
figuration there are N{n— 1) zero modes. Thus there are 
N soft modes whose stiffnesses vanish as we approach the 
fully coUinear configuration. 

For a small displacement away from the coUinear con- 
figuration let 9 be the average of all 9's of all spins. We 
will now argue that for the soft modes the associated sec- 
ond derivative of the potential goes as A ~ 0^, just as for 
the single tetrahedron. 

In a coUinear configuration consider two zero direc- 
tions parameterized by Zq,, that displace the system 
away from the coUinear configuration parallel to a man- 
ifold of cgs. H contains no terms of the form z^, z|, but 
only term z^z^, higher order terms and terms containing 
other displacements. Let us displace the system in the 
(3 direction to say zq/j. Then Za has to stay zero for the 
system to remain in cgs. But around the new point, the 
mode along the a direction would become Zq^z'^ and so 
A ~ Zqp ~ 0^. There are of course all the other modes 
that need to be taken into account but it is reasonable 
that this result will not change. To further check this, 
we expanded the potential around one particular non- 
coUinear cgs, one in which each primitive unit cell has 
the same spin configuration, with the spins displaced by 
an angle 9 from the coUinear configuration. The results 
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are shown in Figure[S^ and[5|3 (with 9 = 0.3). For given 
n and 9 the "spin-wave" bands are the union of b) with 
(n — 2) copies of a). We see that two zero bands, which 
together contain N modes, become nonzero (except for 
special directions that are of measure zero). We further 
find that stiffnesses of these soft modes go as A ~ 0^. 
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FIG. 5: The eigenvalues along chosen directions for a) = 
h) 9 = 0.3 and c) all eigenvalues of the two lowest bands 
ordered by their magnitude 



B. Quantum ground state 

Here we show that the "power-counting" variational 
argument used for the single tetrahedron goes through 
for the full lattice with only slight modifications. We find 
that the ground state wavefunction is localized around a 
collinear configuration for small g and find its spread as 
a function of g and n. 

We consider a variational wavefunction that is local- 
ized near = 0, spreading in angle by 7 in each of the 
~ nN directions available to it. The kinetic energy of 
such a state is ^ gnN/^"^. The typical angle 9 in this 
state is the (Pythagorean) sum of (n — 1) angles of dis- 
placement that are mutually perpendicular to each other 
and each of order 7, so ~ ^^Jn — 1. The kinetic energy 
is ~ Nn?g/9'^. The are N soft modes so the potential 
energy is ~ N9y/g. Minimizing the total energy with re- 
spect to 9 gives the same estimate as before for its typical 
value in the ground state: 



713 ge 



(12) 



Thus we see that the ground state is collinear {9o — > 0) 
in the 5 — > limit. The crossover from collinear ordering 
(small 6*0) to a strongly disordered state can be defined 
as occurring at ^ Ij which puts the crossover in the 
ground state at 



(13) 



C. T>0 

We would like to extend the analysis done for the tetra- 
hedron. What is more complicated here, is that the soft 
modes are not all the same and further they depend in 
some complicated way on the displacement. We found 
that for one particular displacement, their stiffnesses go 
as Xi = ki9^ where fc^'s are shown on Figure Et. To 
proceed we are going to make the following assumption: 
Given some general displacement that is characterized 
some mean 9rms of displacements of all spins, assume 
that this distribution of fc^'s is not going to change very 
much, at least in a qualitative way. 

The probability of finding the system in some cgs, 
treating the perpendicular directions in harmonic ap- 
proximation is given by 



P 



n 



1 



sinh(VA;^) 



where the A^'s are the eigenvalues of the matrix of the 
second derivatives of the potential in the perpendicular 
directions. With our approximations the probability of 
finding the system with some 9 — 9rms is 



P{9)^1[ 



smh{^/Fi9^) 



To find the localization, we need to find the behavior 
of this function. We find that for large n it is sharply 
peaked at 



'V9 



To see this differentiate its logarithm with respect to 9 



d9 



log P{9) 



N{n-2) 



E 



T 



tanh 



for large n. 



and analyze the resulting function for large n. The fc's 
satisfy < ki < 1. For 6*-^ <C 1 the second term is ap- 
proximately N and so the derivative is positive and the 
function is increasing. For 9^ '3> 1 for most of the val- 
ues of i, the term under summation sign is approximately 
\/ki9^ and so the second term is approximately N9^. 
It is a bit smaller because fc^ < 1 but of that order. The 
derivative is zero when 9^ ~ n and is negative for larger 
values of 9 and so as said, the function is peaked at 9f. 
Evaluating the second derivative we find that the func- 
tion is sharply peaked with ratio of the width to 9f being 

^ Therefore the system is localized up to approxi- 
mately 9f just as in the case of the single tetrahedron. 

As n decreases the P becomes more and more dis- 
tributed around zero and at n = 3 the it is only a de- 
creasing function. Nevertheless most of it's weight is still 
in the region approximately < 9f and so the system is 
still localized approximately up to 9f. 



D. Scaling argument IV. NUMERICAL SIMULATIONS 



The previous argument shows that the spins locahze, 
but assumes that the distribution of ki doesn't change 
significantly for other displacements. In this section we 
relax this condition and only assume that A ~ 0^. We 
will show that if the system localizes it depends on g and 
T only through 6*0 = T/y^. 

We want to know the probability that the 9's have rms 
Srms = &r- This is givcu by 

I n ^^(n. sinh( vM{g})/go))-M(yEW - Or) 

The / symbolic, it means to sum over all 

ground state configurations, in the neighborhood of the 
collinear configuration. By changing variables and using 
the scaling property of the eigenvalues it is easy to see 
that 

Thus the P and depends only on 9 /do (and n). Therefore 
the localization depends on g, T only through ^o- 



E. Spin wave theory 



In this section we use imaginary-time path-integral 
quantum Monte Carlo simulations to study the order- 
ing of our quantum rotors on the pyrochlore lattice in 
the Heisenberg case n = 3. This method is based on 
writing the partition function as the trace over all states 
and inserting Nt — 1 additional resolutions of the iden- 
tity to obtain Nt copies ( "time slices" ) of the pyrochlore 
lattice model, coupled ferromagnetically along the imag- 
inary time direction. The case of one time slice is the 
classical model that does not order. We find that an or- 
dered phase does appear already in the case of two time 
slices, although it is restricted to very low temperature. 
As the number of time slices is increased, the ordered 
phase apparently becomes more stable. 

In principle, to get the quantitatively correct behavior 
of our quantum Hamiltonian, we should take the con- 
tinuum (Hamiltonian) limit of A't — > oo (with the ap- 
propriate scalings of Kg and Ki). However, doing this 
extrapolation properly is a computationally demanding 
task that we have not seriously attempted. Instead we 
have simulated primarily the case Nt = 8, mapping out 
a portion of its phase diagram and characterizing the 
phase transition in to the spin-nematic ordered phase. 
We find that the ordering transition is first-order, and 
that the phase diagram is indeed qualitatively similar to 
the crossover diagram we obtained for a single tetrahe- 
dron. We expect that these conclusions are correct also 
in the Hamiltonian limit. 



We now consider a, quadratic, spin wave analysis about 
the collinear cgs to see if a state selection in this approx- 
imation occurs. We note that this should be asymptoti- 
cally accurate at small g and that in the spin system this 
corresponds to the 1/5" computation which does select a 
subset of the collinear states !71. 

A given collinear cgs can be concisely specified, upto 
a global rotation, by writing the spins as Si = rjiZ 
where rji = ±1. Fluctuations about such a configura- 
tion can be parametrized as Si = (xj, ?7i-\/l — x^) where 
Xi = {x^ , . . . ,x"~^). Expanding the square root to sec- 
ond order in a;, the nearest neighbor interaction can be 
written as 



V 



<i.]> 



J nn I 



In every collinear cgs the first term is the same and the 
sum in the last term is ^ rjj — — 2ryi which elim- 

inates all rjs from the last term. Thus the potential 
expanded to the second order around a given collinear 
configuration is completely independent of the choice 
of collinear configuration. Hence any selection beyond 
collinear ordering must come from higher orders in the 
expansion and thus is a weaker effect than one might have 
guessed a priori. 



A. Quantum Monte Carlo 



The system we simulate has partition function 



Z 
S 



Y[dri 



(14) 



Kq Y ^ks ■ nk+i.i - Ki Y, ^k,i ■ rikj (15) 

k,i k,'Ci,j> 

NtT 



NtT 



(16) 
(17) 



where ftk.i is a unit vector at site i in time slice k. The 
nearest neighbor pairs on the pyrochlore lattice are de- 
noted < i,j >, and the relations between the couplings 
here and T and g are shown. The latter are valid in 
the large Nt limit, although we will use them as an ap- 
proximation for finite Nt- The samples are of size 
primitive unit cells (thus 47V^ spins) in each time slice, 
with periodic boundary conditions. 
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B. Order parameter 

The spin-nematic order parameter that we use is thi 
symmetric traceless tensor 

k,i 

where a and (3 run over the three directions in spin space 
Note this order parameter is defined as a sum over al 
of space and imaginary time at one instant during ou 
simulation. To make a dimensionless combination that i; 
sensitive to the ordering, we use the generahzed "Binde 
ratio" 

93 = V6( 3 ) , 

where the average here is over different instantaneous 
measurements of Q during the simulation. This quantity 
is a measure of the degree of collinear spin-nematic or- 
der. It is zero for randomly oriented spins and increases 
as one approaches and enters the ordered phase, tak- 
ing the value one in the well-ordered limit. It is thus a 
good quantity for doing a finite-size scaling analysis of 
the phase transition, as we show below. Note that the 
average value of the third power TrQ^ of the order pa- 
rameter does not vanish in this model. This reflects the 
fact that its Landau theory has a cubic term and thus 
the phase transition is expected to be of first order, just 
as in the case of the isotropic-to-nematic transition in 3D 
liquid crystals. 

C. Phase diagram and the order of the transition 

We first establish the existence of the ordered phase. 
We simulate the system with sizes Nt = 8, = 4, 5, 6, 7 
for Kq = 3 and Ki = 10.1 or Ki = 10.2. In the process, 
for every configuration in the run, we save the "action" 
^0 due to the couplings along the imaginary time direc- 
tion and 5*1 due to the space direction. This allows us 
to obtain estimates of 53 for not only the values of the 
couplings simulated, but also for nearby values of the cou- 
plings by each spin configurations its Boltzmann weight 
^KoSo+KiSi calculating the averages of (73. 

The results for for Kq = 3 and a range of Ki near 
the phase transition are shown in Fig. 5. A crossing of 
the curves for the different sizes is clearly seen at the 
estimated "critical" value of Kic = 10.20. 

Next we find how the order parameter scales with the 
size of the system. We rescale {Ki — Kic) — *■ {Ki — 
Kic)Nx^'^ and tune v to see when the curves of (73 for 
different sizes align. We find that they align well for 
h' = 1/3, as shown in the inset to Figure 5. Thus we 
can write (73 = f{{Ki — Kic)N^) where d = 3 is the 
dimension of the system. This is the scaling expected for 
a first order phase transition. 



1 




0.75' ' ' ' ' ' ' ' ' ' ' 

10.05 10.1 10.15 10.2 10.25 10.3 10.35 10.4 10.45 10.5 10.55 

FIG. 6: The values of ga for sizes A^'^ =4,5,6,7 (as in- 
creases, the curves become steeper). The curves cross at one 
point, which is the estimated location of the phase transi- 
tion. The inset shows qz plotted vs. the scaling variable 
(A'l — Kic)N^ appropriate for a first-order transition. 

To give further evidence of for the first order phase 
transition we calculate the furthest point correlation 
function defined as follows 

2 < {'n-t,x,y,z ■ nt+Nt/2,x+Nj2,y+Ny/2,z+N,/2) > - 

where the average is over all spins at one instant during 
the simulation. If the transition were first order, the finite 
size system at the transition point would jump back and 
forth between the ordered and disordered phase. Thus 
if we plot the histogram of these correlations we should 
see two peaks. These should get sharper as we increase 
the size of the system. The numerical results are shown 
in Figure where we indeed see two peaks, which are 
getting sharper with increasing size. These results, com- 
bined with our expectations from the Landau theory and 
the finite-size scaling of give strong evidence that the 
transition is first order. The place where we are show- 
ing this corresponds to a 5 near the maximum of T^, but 
we expect this first order character will also remain else- 
where on the phase boundary, including at the T — 
quantum phase transition at gc- 

Now that we have explored one point in the phase dia- 
gram, we would like to look at more of them and find the 
shape of the phase transition curve. In order to simplify 
the calculations we will look only for one size = 5, 
Nt — 8 and say that the point is a point of phase transi- 
tion if its qs value is the same as that at the crossing point 
in Fig. 5. Finding a point sufficiently close to the phase 
transition curve we reweight configurations as described 
above to find a point with this value of q^. The results are 
shown in Figure|Sl The range that is readily accessible to 
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Furthest point correlation function 
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FIG. 7: Probability distributions of the furthest-point corre- 
lation function for Nt = 8, Ko = 3 at the values of Ki where 
the two peak heights match. The peaks become sharper and 
the minimum in between deeper with increasing system size, 
as expected for a first-order transition. 



FIG. 8: The estimated phase transition curve for the Nt — S, 
Nx = 5 system. The strong reentrance at low T on high 
g side of the phase diagram is likely an artifact of the very 
coarse imaginary time slicing that is done there with these 
parameters. 



our numerical simulations corresponds to the higher tem- 
perature part of the phase boundary. We see that it has 
quahtatively the same shape as that of our single tetra- 
hedron crossover diagram, and in particular it shows the 
reentrant behavior on varying g at a fixed T (in practice 
this is varying Kq at a fixed Ki). The strong reentrance 
also seen at low T on the high g side of the phase di- 
agram is likely an artifact of the very coarse imaginary 
time slicing that is done there with these parameters. We 
note that the transition temperature is very low, of order 
Tc ~ 1% of J (J = 1 in this paper) and is a manifestation 
of strong frustration in this system. 



V. LARGE N ANALYSIS 

In this section we analyze the system in the limit n 
oo where rotor models are known to be solvable by a 
saddle point method. In the previous analysis we found 
that in this limit the spins don't order at any g and we 
will recover this result. Additionally, we obtain the spin 
correlations in the disordered phase. 

A brief recap of the analysis: We write the parti- 
tion function as a path integral and impose the fixed 
length constraint by introducing a Lagrange multiplier 
field. The resulting action is quadratic in spins which 
can then be integrated out to get an effective action for 
the Lagrange multiplier field with an overall coefficient 
of order n which thus enables a saddle point treatment. 



At T = the saddle point condition takes the form 

where A is the uniform, saddle point value of the La- 
grange multiplier field, a runs over the four bands of the 
pyrochlore lattice and fig are the eigenvalues of the ad- 
jacency (interaction) matrix for a lattice of A'' sites. As 
the the lowest two bands are fiat, with /i'^ independent 
of q, this equation has a solution with A > for any 
nonzero value of g. Hence the system is always quantum 
disordered at n = oo. 

We will not write the resulting correlations of the spins 
explicitly. Instead it is instructive to write the correla- 
tions of the field that is obtained from them as follows. 
Let ft, defined at any site, be a vector pointing from one 
tetrahedron to the other (this can be chosen consistently 
as the tetrahedra surround sites of the bipartite diamond 
lattice). Then define the vector field B°'{x) = n{x)S°'{x). 
The significance of this field lies in the fact that it has 
zero lattice, and thus coarse-grained, divergence in the 
classical ground state manifold. As discussed in Ref. 
this field exhibits dipolar correlations in the T — > limit 
in the classical problem. In our quantum problem, at 
small g and small we find that it exhibits the ground 
state correlations 

{Bt{q)B){^q)) = 

V -f g'^ uo'^ + + gq^ J 
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where we have ignored, for simpUcity, numerical factors 
that appear in front of g. 

These forms reduce at equal times and 5 = to the 
dipolar forms derived in At nonzero g we see that 
the correlations decay exponentially in space with a cor- 
relations length of order g^^l"^ and exponentially in time 
with a gap of order g. Observe that the n = 00 problem 
does not contain any trace of the nematic ordering that 
exists at finite n. Individual terms in the 1/n expansion 
for the correlation function can be seen to be well be- 
haved and so it will take an analysis of the series to 
reproduce the instability that we have obtained earlier. 
This is a challenge for future work. 

VI. CONCLUSIONS 

We have fully explored the local ordering of a single 
tetrahedral unit of four neighboring rotors. They or- 
der collinearly, with two rotors pointing along one direc- 
tion and the other two in the opposite direction. Thus 
an axis of ordering is chosen, and since in the full py- 
rochlore lattice this tetrahedron shares one rotor with 



each of its neighbor tetrahedra, this axis of ordering 
is uniform throughout the lattice in the ordered phase, 
which thus has spin-nematic, or collinear order. We have 
demonstrated this spin-nematic ordering within a quan- 
tum Monte Carlo simulation. There might also be long- 
range order in which rotors point in which direction along 
this axis, and/or in which pairs of rotors are parallel or 
antiparallel. We have not yet been able to determine 
whether our ordered phase has any of these latter types 
of sublattice order in addition to its collinear order. We 
find that the collinear order shows up quite robustly in 
our quantum Monte Carlo simulations, while any other 
long-range correlations that might be present appear to 
be much weaker and difficult to detect, if they are indeed 
there. 

We note that the considerations outlined here apply 
straightforwardly to the the two-dimensional checker- 
board lattice which is the planar analog of the pyrochlore 
lattice. By contrast an analysis of the rotor models on 
the kagome lattice will require a fresh analysis. 

We thank the NSF for support through MRSEC grant 
DMR-0213706. 
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